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Abstract. The true nature of the progenitor to GRBs remains elusive; one characteristic that would constrain 
our understanding of the GRB mechanism considerably is gamma-ray polarimetry measurements of the initial 
burst flux. We present a method that interprets the prompt GRB flux as it Compton scatters off the Earth's 
atmosphere, based on detailed modelling of both the Earth's atmosphere and the orbiting detectors. The BATSE 
mission aboard the CGRO monitored the whole sky in the 20 keV - 1 MeV energy band continuously from April 
1991 until June 2000. We present the BATSE Albedo Polarimetry System (BAPS), and show that GRB 930131 and 
GRB 960924 provide evidence of polarisation in their prompt flux that is consistent with degrees of polarisation 
of II > 35% and n > 50% respectively. While the evidence of polarisation is strong, the method is unable to 
strongly constrain the degree of polarisation beyond a systematics based estimation. Hence the implications on 
GRB theory are unclear, and further measurements essential. 
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1. Introduction 

Whatever the mechanism, polarisation in the prompt 
gamma-ray flux from GRBs is evidence of strong mag- 
netic fields within the burst. Theories on the GRB pro- 
duction mechanism can be constrained by different de- 
grees of linear polarisation. For large degrees of po- 
larisation, IT « 80%, either shock accelerated syn- 
chrotron emi ssion or a tuned Compton-drag model is the 
most likely llLazzati et all l2004t ICoburn fc Boggj 120031 
iLvutikov et al.ll2003|) . For intermediate degrees of polari- 
sation (20%< LT < 60%) two electromagnetic models have 
emerged that involve either synchrotron emission as the 
dominant source of radiation or as the result of viewing: 
the burst from just outside the edge of the jet l|Granotl 
l200fl iGhisellini fc Lazzatlll99^ . Low degrees of polari- 
sation can be a result of either a hydrodynamic model or 
from flux with a higher degree of polarisation experiencing 
partial depolarisation. 
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The high degree of linear polarisation initially reported 
in the prompt flux of GRB 021206 [deg ree of polarisation, 
n = 80±20 %.ICoburn fc BogpfeOOaV i with the RHESSI 
experiment l|McConnell et alJl2002l) has stimulated much 
interest in the implications this has on GRB theory. 
Though GRB astronomy has the advantage of large fluxes, 
any polarimetric measurements are still dominated by sys- 
tematic effects that can only be properly quantified by 
careful modelling. The importance of correctly evaluating 
the systematic effects is paramount in any measurement 
of GR B polarisation as was emphasised bv IWigger et alJ 
(2004) in the thorough re-anal y sis of the initi a l RH ESSI 
result l|Coburn fc Boggj |2003|) . IWigger et all (|2004 re- 
analysed the degree of polarisation to be II = 41+44%. 
This implies that many of the production mechanism the- 
ories are now similarily competitive. 

Besides the RHESSI instrument l)Lin et all 12002^ . 
there are several other instruments that could be suit- 
able gamma-r ay polarimete r s: IB IS and SPI aboard 
INTEGRAL jWinkler et alJ 12003^ have the as yet 
untested capabili ty of detecting gamm a-ray polarisation; 
COMPTEL l|Schoenfelder et alJll993ft and BATSE, both 
aboard the Compton Gamma-Ray Observatory (CGRO) 
Irishman et alJfl989gJ) are also possible candidate instru- 
ments. For Compton polarimetry, the instrument is re- 
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quired to scatter the photons and detect the modulation 
induced by the polarisation-sensitive different ial cross- 
sectio n for Compton scattering (equation^), see ijLei et alJ 
119971) for a review of Compton polarimetry in gamma-ray 
astronomy. 
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where r 2 , is the classical radius of the electron, Eq the 
energy of the incident photon, E' the energy of the scat- 
tered photon, 9 the scatter angle and (f> the azimuthal 
angle relative to the polarisation direction. A problem in 
both COMPTEL and RHESSI is that each instrument 
has to both scatter and detect the linearly polarised pho- 
tons. The scattering effective area, A e ff, is small, result- 
ing in a reduction in the number of double events re- 
quired to detect the polarisation-induced azimuthal mod- 
ulation. COMPTEL also had the problem that to act as 
a polarimeter th e double-event rate had to be quite high 
ijLei et al.lll996|) . and at these high event rates vital kine- 
matic information was lost in the telemetry. This is also 
expected to affect the polarimetry mode in IBIS to some 
extent. 

To both scatter and detect the polarisation-induced az- 
imuthal modulation, an instrument has to be either highly 
event selective and tracking sensitive or requires a large 
scattered flux, larger than the maj ority of GRBs can sup - 
ply to RHESSI and INTEGRAL (jSwinvard et al.lll99l|) . 
The problem of tracking and selecting the events most effi- 
ciently is currently being addressed in missions such as the 
Tracking and Imagin g Gamma-Ray Experiment ( TIGRE) 
I O'Neill et alJ Il99fil) an d the Gamma-RAy Po larimeter 
Experiment (GRAPE) l)McConnell et alJl2003|) . BATSE 
however is in the unique position to address the prob- 
lem of increa sing A e f f for scattering t he incident gamma- 
ray photons. iMcConnell et al.l <|l994^ proposed that po- 
larimetry was possible by combining the large A e ff of 
the Earth's atmosphere for Compton scattering with the 
4-7T steradian field-of-view of the eight BATSE modules. 
Unfortunately the response of the CGRO to non-planar 
flux was difficult to handle analytically due to the un- 
collimated nature of BATSE and so could only be prop- 
erly evaluated with modelling techniques. The modelling 
technique, called Mass Modelling, has been developed 
at South ampton and has b een employed successful ly on 
BATSE JWestmord (l2002h : IShaw et all J2003L l2004h and 
references therein) IN TEGRAL l|Ferguson et alJ l2003|) 
and Swift ijWillisl l2002j) . The complex analytical situa- 
tion and systematic effects resulting from immersing the 
CGRO in the GRB flux scattered off the Earth's atmo- 
sphere can be modelled using this Mass Modelling ap- 
proach, enabling polarimetry measurements with BATSE 
to take full advantage of the large A e f f that the Earth's 
atmosphere provides. 

Here we present the technique of alb edo Compton po- 
larim etry from the Earth's Atmosphere I McConnell et alJ 
Il994h combined with the BATSE Mass Model (BAMM) 



l|Shaw et alJ l2003|) . the BATSE Albedo Polarimetry 
System (BAPS). 

2. The BATSE Albedo Polarimetry System 
(BAPS) 

The BATSE mission aboard CGRO monitored the 
whole sky in the 20 keV - 1 MeV energy band continu- 
ously from April 1991 until June 2000. I t was designed 
primar ily to detect and locate GRBs i Fishman et alJ 
Il989b|) . BATSE consisted of eight uncollimated 2025cm 2 
Nal(Tl) Large Area Detectors (LADs), which had a 
combined field of view of 4tt steradians. Also positioned 
below each LAD was a Spectral Detector (SD) made 
of 127 cm 2 Nal(Tl). The SDs were optimized for en- 
ergy coverage and resolution, and were operated over 
the range lOkeV - 100 MeV. An illustration of the 
CGRO and a BATSE module can be seen in Figure 
n BATSE was an extremely successful experiment 
and the current burst catalogue has over 8000 triggers 
(http://cossc.gsfc.nasa.gov/batse/BATSE_Ctlg/index.html). 
The discovery that GRBs are located at cosmological 
distances is largely due to data collected with BATSE 
l|Meegan et al.lll992|) . 

BAPS uses both the BATSE Mass Model (BAMM) 
and a model of the Earth's atmosphere. The system is 
highly CPU intensive and combining both models into 
one model is impractical due to the relative scales. The 
model of the Earth's atmosphere was written in GEANT4 
l|Ravndal et alJll993|) . in order to include the polarised 
Compton scattering process. It consists of a ball of ra- 
dius 6503 km [Earth radius + 125 km] made of an inte- 
grated atmospheric composition of constant density [p = 
8.33 x 10~ 5 g.cm- 3 , 38% 2 : 62% N 2 ], the effects of this 
assumption are discussed in section E.4.31 The BATSE vol- 
ume was an imaginary sphere of radius 200 km at an al- 
titude of 500 km. Figure [3 shows the scale of the Earth 
albedo model alon gside a schematic of BAMM. BAMM 
l|Shaw et alJl2003j) was written in GEANT3 to dynami- 
cally flat-field the var iable background for the purposes of 
occultation imaging l(Harmon et alJl2002l) and the subse - 
quent production of an all-sky survey IjShaw et alJ 120041) . 
BAMM reproduces the instrument performance to a level 
of accuracy that is impossible through analytical tech- 
niques alone. For a full review of the mass model ling ap- 
proach and its applications, see lDean"et al. (20o3). 

BAPS simulates the GRB flux as it scatters off the 
atmosphere; the distribution of this scattered flux is 
recorded as it passes through a volume equivalent to where 
the CGRO was at the time of the burst. Polarised flux 
preferentially scatters perpendicular to the direction of 
the polarisation vector, so any distribution produced as 
a result of polarised flux will appear as an anti-phased 
excess towards the limbs of the Earth, when viewed from 
the direction of the burst. The recorded distribution scat- 
tered from the atmosphere is then used as an input to 
the BAMM Monte Carlo and the count-rates simulated 
for each LAD. This process is carried out for an unpo- 
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larised burst and for a fully polarised burst, for a range of 
polarisation angles. The simulated LAD count-rates are 
then compared to the burst count-rates in the BATSE 
data. The fitting procedure varies the percentage polari- 
sation and the polarisation angle. The Z statistic, defined 
by Z = — R i- R 2 where R and t are the count-rates and 

observation times respectively, is used to compare the two 
count-rates: simulated and recorded. 

BATSE uses two formats of data, CONTinuous 
(CONT) and DISCriminator Large Area (DISCLA). 
DISCLA data contains 4 broad energy channels with a 
binning width of 1.024 s. CONT data contains 16 energy 
channels with a binning width of 2.048 s. Due to BAPS 
only being sensitive to the Compton domain and below, 
CONT data is used to optimise this sensitivity. The sac- 
rifice of better timing resolution is not a problem until 
time-evolving polarimetry is considered. 



3. Burst Selection: GRB 930131 and GRB 960924 

The selection criteria for a suitable burst are that the 
burst has to be short, strong, have a geocentre to burst 
angle, 77, as close to 180° as possible (i.e. to have BATSE 
directly between the Earth and the burst) and for it to 
have no obvious evolving features. This selection reduces 
any background variations, maximise the polarisation- 
induced modulations and avoid the possibility of repro- 
cessing at the source reducing the degree of polarisation, 
perhaps due to reverse shocks and scattering at the source. 
Suitable bursts were identified via a plot of burst strength 
against burst to geocentre angle, rj (Figure |3J). Two of 
the strongest and most suitably placed BATSE triggers 
were identified as being suitable for this study: trigger 
2151 (GRB 930131) and trigger 5614 (GRB 960924). The 
remaining triggers were either too long or had clear evolv- 
ing features in the light curve, and so were discarded at 



Table 1. A table of the relevant parameters to BAPS. The 
spectral parameters used were obtained usi ng a standard 
spectral characterisation ijBand et alJll993]) . 



GRB 


930131 


960924 


Trigger 


2151 


5614 


O 


170.7 


162.4 


CONT bins (2.048 seconds) 


2 


4 


Band model a 


-1.2 


-0.3 


Band model f3 


-2.5 


-9.9 


Band model Et rea k (keV) 


439 


251 



this stage. The parameters of these two selected bursts are 
displayed in tabled 




Fig. 3. All BATSE triggers plotted to identify the six most 
suitable candidates for polarimetry. Triggers 2151 and 
5614 were initially selected with the remainder discarded 
due to their particularly long and structured lightcurves. 
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Fig. 2. (a) This shows the relative scale of the Earth, its atmosphere (up to 125 km) and the virtual BATSE volume 
for the GEANT4 polarisation sensitive simulations. (b) BAMM contains detailed structural, geometrical and material 
infomation for BATSE, EGRET, COMPTEL, OSSE and the CGRO spacecraft. Cut-aways of the BATSE module 
(top right) and the whole CGRO (bottom right) show the internal detail. To illustrate the crucial components that 
will deflect or absorb gamma-ray photons, some external volumes have been made transparent. 



4. Analysis 

4.1. Simulations 

The first stage of BAPS is to produce the input distribu- 
tions for the Mass Model. Fully polarised flux from the 
GRBs (with polarisation angle incremented every 15°), 
and an unpolarised burst flux were simulated incident on 
the atmosphere. 0.03% of the incident flux was seen to 
reflect and pass through the imaginary bubble containing 
the CGRO. The albedo angular distribution for polarised 
and unpolarised flux for both the GRBs can be seen in 
Figure 

BAMM was then used to simulate the effect these dis- 
tributions of albedo flux have on the count-rates in the 
LADs during each burst. The direct GRB flux is also 
a simulated component. As seen in the albedo spectra 
(Figure EJ) this technique is only sensitive in the Compton 
domain and so the energy range of the fitting is restricted 
to LAD CONT channels 3-6 (31-98 keV). The Compton 
domain extends up to higher energies but these restric- 
tions equated to 67% and 62% of the albedo flux for 
GRB 930131 and GRB 960924 respectively. 

4.2. LAD selection and fitting 

The LADs selected are divided into two types, the first 
are used to normalise the model count-rates for the direct 
and albedo fluxes and the second are LADs pointed such 
that the albedo count-rate varies considerably depending 
on the degree and angle of polarisation. The remaining 
LADs are ignored if they have low simulated count-rates 
and low levels of polarisation-induced modulation, thus 
reducing the systematic and counting errors in the system 
and avoiding adding unconstrainable parameters into the 



Fully Polarised Flux Unpolarised Flux 




Fully Polarised Flux Unpolarised Flux 




Fig. 4. The simulated angular distribution of gamma- 
ray flux reflected off the Earth's atmosphere during 
GRB 930131 and GRB 960924 for 100% polarised and un- 
polarised fluxes. Contours are the 0, 5, 10, 15, 20 and 
25 mean counts per pixel (4.5° x 4.5°) for 5xl0 7 incident 
GRB photons. This pixelated distribution is used as the 
input to BAMM. 



fitting process. The component direct from the burst is 
normalised using the count-rate from the most burstward 
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Fig. 5. The spectra of the albedo component of 
GRB 930131 and GRB 960924. The dotted lines represent 
the domain in which the data is compared: 31-98 keV 
(LAD CONT channels 3-6). These limits equate to 67% 
and 62% of the albedo flux respectively. 

LAD. The simulated count-rates are fitted to the recorded 
count-rates with the normalisation of the direct flux fixed 
and the albedo flux, degree and angle of polarisation as 
free parameters. The ambient background is removed by 
taking a mean count-rate from before and after the burst. 
The fit statistic used is a simple Z-test for comparing the 
simulated and recorded count-rates. 

The criteria for selecting the LADs is best demon- 
strated by looking at the simulated count-rate components 
for each LAD, see tables and 03 The figure-of- merit con- 
sists of the percentage of simulated albedo flux, A, in each 
LAD compared to the total flux, the maximum variation 
due to GRB polarisation angle and the level of albedo flux 
compared to the LAD with the maximum albedo flux. For 
each burst one LAD is selected to normalise the direct 
burst flux, one to assist in the normalisation of the albedo 
flux and two to measure polarisation-induced modulation. 
For different bursts the number of LADs selected may vary 
according to their figure-of-merit. 

4.2.1. GRB 930131 

The LADs selected are LAD1, LAD2, LADS, and LADl 
LAD1 and LAD2 have the largest albedo count-rate and 
are most sensitive to polarisation. LAD4 is used to set 
the normalisation for the direct flux and LAD3 is used to 
assist the normalisation of the albedo flux. 

4.2.2. GRB 960924 

The LADs selected are LAD0, LAD2, LADl and LAD6. 
LAD 6 is used to set the normalisation for the direct 



Table 2. GRB 930131 LAD figure-of-merit. The LADs 
used in the fitting are marked in bold. 



LAD 



A + Direct 



% Var 



% 



% 



LAD0 
LADl 
LAD2 
LAD3 
LAD4 
LAD5 
LAD 6 
LAD7 



1 

49 
70 
73 
<1 
16 
<1 



45 
38 
17 
7 

58 
25 
38 
31 



38 
59 
60 
100 

1 

2 

9 
17 



Table 3. GRB 960924 LAD figure-of-merit. The LADs 
used in the fitting proceedure are marked in bold. 



LAD 



A + Direct 



% Var 



% 



% 



LADO 

LADl 
LAD2 
LAD3 
LAD4 
LAD5 
LAD6 
LAD7 



31 
51 

8 

7 

3 

4 

<1 
4 



11 
10 
30 
35 
36 
44 
21 
41 



100 

29 
57 
14 
62 
13 
8 

<1 



flux. For this orientation of the CGRO, the polarisation- 
induced variation appears in all LADs. LADO is used 
to primarily assist the normalisation of the albedo flux 
whereas LAD2 and LAD4 show the largest albedo flux 
and have > 30% variability due to polarisation angle and 
so are used to test any polarisation-induced modulation. 

4.3. Results 

The selected LADs are simultaneously fitted to the CONT 
count-rates and the best fit for each incident polarisation 
angle is shown in Figures and [H] The significance of 
each of these polarisation angle fits can be taken from 
Figures and 13 the values for which were taken from a Z- 
test comparing the total count-rates in the LADs from the 
simulation and CONT data using the parameters supplied 
by the best fit for each polarisation angle. Also included 
in these plots of the Z statistic are the corresponding II 
values. 

GRB 930131 shows evidence of a high degree of polar- 
isation. In this case, using simulated data with the polari- 
sation angle at 90 degrees to the actual polarisation angle, 
the best fit will be to the unpolarised simulated flux, due 
to the phased nature of the excess at the Earth's limb. 
Hence for this burst the degree of polarisation will only 
move towards a higher value when the fit is improved by 
the polarised GRB. Though this is seen in Figure the 
sharp change in the degree of polarisation compared to 
the smoother change in the Z statistic suggests that this 
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Fig. 6. A comparison of the fitted simulated count-rates with the LAD CONT data for GRB 930131. The CONT data 
is shown as a dashed line, the error on that count is shown as solid lines. The best fit for each polarisation angle is 
shown with the relevant counting error from the Monte Carlo. LADS and LAD4 simulations are expected to match 
the recorded data, independent of polarisation. LAD1 and LAD2 are sensitive to polarisation in the GRB's prompt 
flux and will show an improvement in fit for angles that corresponds to the burst polarisation angle. 



100% 




Polarisation Angle (degs) 

Fig. 7. For GRB 930131, the Z statistic (solid line) for 
each of the best fits from Figure The statistic becomes 
significant (95% confidence) between polarisation angles of 
163° and 9° and this corresponds to a high (>90%) degree 
of polarisation (dashed line). The Earth's North/South 
equivalent angle is at 5°, see section 14. 4. II 
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0% 

Polarisation Angle (degs) 

Fig. 9. For GRB 960924, the Z statistic (solid line) for 
each of the best fits from Figure El The statistic becomes 
significant (95% confidence) between polarisation angles of 
58° and 64° and this corresponds to a high (100%) degree 
of polarisation (dashed line). The Earth's North/South 
equivalent angle is at 47°, see section 14.4.11 



technique is not particularily sensitive to the degree of po- 
larisation and can initially only be used as evidence of po- 
larisation in the prompt GRB flux and that GRB 930131 



is best fit by a fully polarised model (though see section 
14.51 for an estimate of II). 



D.R. Willis et al.: The BATSE Albedo Polarimetry System (BAPS) 



7 




Polarisation Angle (degs) 



Polarisation Angle (degs) 




Polarisation Angle (degs) 



LAD 6 










j j cjoNljata j 








l j j ~ a 







Polarisation Angle (degs) 



Fig. 8. A comparison of the fitted simulated count-rates with the LAD CONT data for GRB 960924. The CONT 
data is shown as a dashed line, the error on that count is shown as solid lines. The best fit for each polarisation 
angle is shown with the relevant counting error from the Monte Carlo. LAD6 simulations are expected to match the 
recorded data, independent of polarisation. LAD2 and LAD4 are sensitive to polarisation in the GRB's prompt flux 
and will show an improvement in fit for the specific angle that the burst is polarised to. LADO is largely insensitive 
to polarisation and so should be consistent with the CONT data regardless of polarisation. 



GRB 960924 also shows evidence of a high degree of 
polarisation. The variation in the Z statistic in this case 
shows a clear dip between 35° and 120° and is flat for the 
remainder of angles. The sharp increase in the degree of 
polarisation compared to the dip in the Z statistic clearly 
shows the angles for which a polarised burst fits the CONT 
data better than an unpolarised burst. Again this is evi- 
dence of polarisation, that is best fit by a fully polarised 
model. 

4.4. Systematic Effects 

These results are expected to contain a series of system- 
atic effects and a discussion and evalution of these effects 
follows. The LAD selection by way of a figure-of-merit for 
each burst is expected to minimise these. 

4.4.1. North-South Albedo Gamma-Ray Excess 

The atmosphere itself is a strong source of gamma radi- 
ation. Cosmic ray interactions in the atmosphere radiate 
gamma-rays isotropically. This emissivity depends solely 
on the intensity of the incident cosmic rays and the at- 



mospheric depth IjDean et al.lll989l) . An angle- averaged 
spectr um of this albedo flux is also presented bv iGehrelsl 
This albedo flux will be concentrated at the 
Earth's limb as the largest portion of the atmosphere is 
being observed. Along the limb, the strongest source of 
albedo flux will be towards the polar regions, as the rigid- 
ity of the ge omagnetic field defle cting the cosmic rays is at 
its weakest ((Smart fc Sheall994h . This will result in an ex- 
cess that will be reminescent of the excess that BAPS will 
be looking for. During these bursts BATSE was at an incli- 
nation of -16° and -10° for GRB 930131 and GRB 960924 
respectively. The emission from the south pole would have 
been visible to BATSE in both cases. The cosmic ray in- 
duced albedo flux would therefore falsely suggest a polari- 
sation angle of 5° and 47° respectively, which corresponds 
to the inferred polarisation angle for GRB 930131 but not 
for GRB 960924. However, GRB 930131 should not be dis- 
counted as a possible source of prompt polarised flux as 
any excess produced from the South pole should be re- 
moved during the background subtraction within the fit- 
ting proceedure. 
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4.4.2. East-West Cosmic Ray Anisotropy 

Another effect of the geomagnetic field is the cosmic ray 
east-west anisotropy, from the direct detection of cosmic 
ray particles. The rigidity of the geomagnetic field indi- 
cates the ability of the field to fully deflect a charged 
particle. Faraday's Law defines the direction that a cos- 
mic ray is deflected in, and results in an excess in the 
westerly direction. This would correspond, in this case, to 
a falsely inferred polarisation angle of 95° and 137° for 
GRB 930131 and GRB 960924 respectively. These polari- 
sation angles do not show an improvement in the fit and 
so this anisotropy is not responsible for a positive detec- 
tion. Similar to the polar albedo excess discussed in sec- 
tion any anisotropy should not be seen due to the 
background removal and it is not expected to change on 
the order of seconds. Moreover, any level of background 
due to cosmic ray interactions will be low as each BATSE 
LAD had a plastic scintillator vetoing cosmic ray events. 
Finally, it is estimated that both the North-South albedo 
excess and the East- West cosmic ray anisotropy are ex- 
pected to be extended by 20-30°, leaving few of the poten- 
tial polarisation angles free from these systematic effects. 
For weaker bursts, these North-South and East- West re- 
gions may give false detections. 



(38% O2 : 62% N 2 ,). This is an approximation that will 
introduce a systematic error but aid in the simplicity and 
speed of BAPS. The constant density was acheived by av- 
eraging over the first 125 km of the atmosphere. To eval- 
uate the effectiveness of this approximation, a single col- 
umn element of this atmospheric model was compared to 
a graded (in density and composition) 12 layered model of 
the atmosphere. The rate of back-scattered events (scatter 
angle >90°) and resultant direction of the photon emerg- 
ing from the atmosphere were recorded for both models, 
as seen in Figure EH For the layered model, the major- 
ity of the scattering occurs in the 20-50 km region, and 
the 0-125km integration reproduces 99.7% of these back- 
scattered events. The back-scattered zenith angle distri- 
bution for an on-axis beam of photons and one at 40° 
off-axis was also compared for the two models. There is a 
slight excess around the scatter angle of 160° and below 
for the approximated model. However, it is clear that this 
discrepancy between the two models will not effect the 
BAPS fitting when considering the uncollimated nature 
of the BATSE LADs. Furthermore, the 0.3% change in 
the albedo count-rate will be accounted for in the fitting 
routine. The diametric excess that polarisation will induce 
will not be effected. 



4.4.3. Modelling Software 

Combining two GEANT models will introduce a level of 
systematic error but due to the relative scale of the two 
models this is unavoidable. The GEANT 3 model of the 
CGRO is proven to work to an acceptable level of ac- 
curacy l|Shaw et al.ll2003h . However, the GEANT 4.6.1 
model of the Earth's atmosphere may introduce some 
problems. GEANT 4 is a new suite of software and is still 
under development, requiring verification at each stage. 
One area that requ ires modification was highlighted by 
iMizuno et alJ l(2004h : the modulation produced from scat- 
tering a polarised beam. For their GEANT model the un- 
modified version of GEANT4 was seen to scatter a Crab 
(25-200 keV) spectrum with around 30% variation with 
scattering angle, whereas the beam-test verified and mod- 
ified version of GEANT was seen to give a modulation 
of around 64%. The models are entirely different but this 
suggests that the level of modulation in each LAD due 
to polarisation in BAPS could be up to half of what 
it should physically be. This is undoubtably a system- 
atic error which will degrade the sensitivity of BAPS but 
not one that would give a false indication of polarisation. 
Future runs will inc lude the modifications suggested by 
IMizuno et"aH l)2004|) and it is expected that BAPS will 
become far more sensitive as a result, and will give a bet- 
ter indication as to the degree of polarisation. 

4.4.4. Constant Density Atmospheric Model 

The model of the Earth's Atmosphere consists of a ball 
of constant density (8.33 x 10 3 g.cm -3 ) and composition 
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Fig. 10. A comparision between using a 12 layered graded 
atmospheric model and a constant density and compo- 
sition integrated model. The incident photon energy is 
200 keV. The upper plot shows the back-scatter rate with 
altitude and the lower two histograms show the albedo 
zenith distribution. 
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This depth-dependent scattering rate also discards the 
need for including the Earth itself into the model. A pho- 
ton interacting such that it scatters in the atmosphere to- 
wards the direction of the CGRO but must travel through 
the Earth to reach the CGRO volume will certainly be ab- 
sorbed by the atmosphere in the model, and so discarded. 

4.5. Estimating the Degree of Polarisation 

The physical implications of polarised prompt flux from 
GRBs hinges on II, the percentage of the flux that is po- 
larised. As BAPs normalises the levels of albedo flux and 
involves comparisons of count-rates, the degree to which 
these GRBs are polarised cannot be measured directly. 
However, an estimate of II can be achieved through re- 
stricting II to a maximum value and inspecting the change 
in the Z fitting statistic. The data is refitted with restric- 
tions applied to II, as shown in Figures ITT1 and IT21 
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Fig. 11. The CONT data for GRB 930131 is compared to 
the simulations with II restricted to a range of maximum 
values between 0% and 100%. The confidence levels in the 
Z statistic are shown. 



Though the significance of the fit for GRB 930131 de- 
creases as the maximum allowed value for II decreases, 
the fit is still significant down to a II values of around 
70%-80%. However, the fit for GRB 960924 suffers more 
severely with any change to n and the significance of the 
result is lost at a level of 90% polarisation. This value for 
n is unphysical and so the systematics discussed in section 
14.4.31 must be considered. As it is quite possible that the 
polarisation-induced modulation is underestimated by up 
to a factor of 2, the value of II may well be overestimated 
by a similar figure. Therefore the degree of polarisation 
in GRB 930131 is consistent with II ranging between 35% 
and 100%; and GRB 960624 is consistent with II ranging 
between 50% and 100%. 
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Fig. 12. The CONT data for GRB 960924 is compared to 
the simulations with II restricted to a range of maximum 
values between 0% and 100%. The confidence levels in the 
Z statistic are shown. 

5. Conclusions 

BAPS has been presented here in the context of 
GRB 930131 and GRB 960924, both of which show a high 
degree of polarisation. The caveats to this result have been 
discussed and no systematic explanation can be found for 
these positive detections of polarisation. The GRB selec- 
tion procedure is initially highly restrictive, however there 
are several bursts that are also suitable to be analysed by 
BAPS in future studies. The initial stage of the BAPS 
analysis (modelling of the albedo flux) was run for a va- 
riety of angles taken from figure and it was seen that 
the difference between the polarised and unpolarised cases 
was statistically significant up to 77 < 120°, suggesting 
that BAPS is primarily restricted in its sensitivity by the 
number of LADs that are illuminated by the GRB albedo 
flux and not the off- axis angle of the burst. This increases 
the number of GRBs that can be analysed with BAPS, 
however since the senstivity is governed by the number of 
illuminated LADs, it removes the possibilty of analysing 
the systematics of a suitably off-axis GRB that is unable 
to exhibit any polarisation effects. 

There are areas in BAPS that would benefit from fur- 
ther augmentation: modification to the current GEANT4 
release, an analytical model of the scattering off the 
Earth's atmosphere and a full analysis of the systematic 
errors to enable a quantitative estimate to the degree of 
polarisation and enable the inclusion of more LADs into 
the fitting proceedure. 

Though the degree of polarisation is an estimate from 
these results, there is distinct evidence that polarisation 
exists in prompt GRB flux and to a degree that is large 
enough to be detected by viewing the Earth's atmosphere 
with the BATSE detectors. The implications of this on the 
GRB mechanism are that a large scale magnetic field is 
present in the emitting region, conclusions beyond this will 
remain uncertain until future measurements can be made. 
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Viewing the prompt flux is crucial to future polarimetry 
measurements, something that is only currently possible 
through a large scale survey of the BATSE catalogue with 
BAPS or through an on-axis INTEGRAL burst such as 
the recent strong burst GRB 041219. 
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